Have Superkamiokande Really Measured the Direction of the Atmospheric Neutrinos 
which Produce Fully Contained Events and Partially Contained Events ? 
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Quasi Elastic Scattering (QEL) is the dominant source for producing both Fully Contained Events 
and Partially Contained Events in the Superkamiokande(SK) detector for the atmospheric neutrinos, 
in the range ~0.1 GeV to ~10 GeV. In the analysis of SK events, it is assumed that the direction 
of the incident neutrino is the same as that of the detected charged lepton. In the present letter, 
we derive the distribution function for the scattering angle of the charged leptons, their averaged 
scattering angle and their standard deviation due to QEL. Then, it is shown that the SK assumption 
for the scattering angle of the charged leptons in the QEL is not valid. Further, we examine the 
influence of the azimuthal angle of the charged leptons over their zenith angle. As the result, 
we conclude that the zenith angle distribution of the neutrino under the SK assumption does not 
reflect the real zenith angle distribution of the atmospheric neutrino which produces Fully Contained 
Events and Partially Contained Events. This result has clear implication for attempts to detect 
neutrino oscillations from the analyses of Fully Contained Events and Partially Contained Events 
in Superkamiokande. 
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I. INTRODUCTION 

In the experiment in which they try to detect the 
neutrino oscillation, by using the size of the Earth 
and measuring the zenith angle distribution of the 
atmospheric neutrino events, such as, Superkamiokande 
experiment [1] - hereafter, simply SK it is demanded 
that the measurements of the direction of the incident 
neutrino are being carried out as reliably as possible. 
Among the experiments concerned on the neutrino 
oscillation, the analysis of Fully Contained Events in 
SK is regarded as mostly ambiguity-free one, because 
the essential information to extract clear conclusion is 
stored inside the detector. In SK, they assume that 
the direction of the neutrino concerned is the same as 
that of the produced charged lepton (hereafter, simply 
SK assumption) [2,3]. However, the SK assumption does 
not hold in the just energies concerned for neutrino 
events produced inside the detector, which is shown later. 



and Parially Contained Events (single ring events) are 
analysed, quasi elastic scattering of neutrino interac- 
tion(QEL) is the dominant source for the atmospheric 
neutrino concerned [4] 



II. THE DIFFERENTIAL CROSS SECTION FOR 
QEL AND THE SCATTERING ANGLE OF THE 
CHARGED LEPTON 

The differential cross section for QEL is given as 
follows [5]. 
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In the energy region where Fully Contained Events The signs + and - refer to ^Ve) and f^( e ) for charged 
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FIG. 1: Distribution function for the scattering angle 9 S of 
the muon for muon-neutrino. 

current(c.c) interaction, respectively. The Q 2 denotes 
four momentum transfer between the incident neutrino 
and the charged lepton. As for details of other symbols, 
see the text [5]. The relation among Q 2 and E v , the 
incident energy of neutrino, Eg, the energy of the emitted 
charged lepton ((anti)muon or (anti)electron) and 9 S , the 
scattering angle of the charged lepton, is given as 

Q 2 = 2E v E e (l-cos6 s ). (2) 
Also, the energy of the charged lepton is given by 
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For a given energy E v of the incident neutrino, we 
randomly sample Q 2 through the Monte Carlo procedure 
from Eq. (1). Subsequently we obtain the scattering 
angle 8 S of the charged lepton concerned by Eqs. (2) 
and (3). Thus, we obtain the distribution functions for 
scattering angle of the charged lepton. In Fig. 1, we give 
such distribution functions for different incident neutrino 
energies. Through such a Monte Carlo procedure, we 
obtain the average scattering angles and their standard 
deviations, too. We give them in Table 1. It is shown 
clearly from the figure and the table that the average 
scattering angles largely deviate from the direction of 
the incident neutrino, being accompanied by rather large 
standard deviations and consequenly we cannot neglect 
the scattering angle in the energy region where SK was 
interested in, say ~0.1 GeV to ~10 GeV. 



TABLE I: The average values < 9 S > for scattering angle of 
the emitted charged leptons and their standard deviations a B 
for various primary neutrino energies E v 
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angles of the charged particles in QEL play a decisive 
role in the determination of their zenith angles as well as 
the translation from Fully Contained Events to Partially 
Contained Events (vice versa) which are mentioned later. 

In order to examine the influence of the azimuthal an- 
gle of the charged leptons over their zenith angle, let 
us denote the direction cosines of the incident neutrino 
(£, to, n) and denote the scattering angle of the charged 
lepton, 9 S , and the azimuthal angle, <fi, with regard to 
the axis of the incident neutrino. Then, (£ r ,TO r ,n r ) , the 
direction cosines of the charged lepton which correspond 
to (£, m, n) are given as 
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IV. MONTE CARLO PROCEDURE, TAKING 
ACCOUNT OF THE EFFECT OF THE 
AZIMUTHAL ANGLE AND DISCUSSIONS 



III. THE INFLUENCE OF THE AZIMUTHAL 
ANGLE OF THE CHARGED LEPTON OVER 
THEIR ZENITH ANGLE 

In addition to the scattering angle of the charged 
leptons, it should be cmphaized that the azimuthal 



By using Eq. (4), we carry out a Monte Carlo 
calculation to examine the influence of the azimuthal 
angle of the charged leptons over their zenith angle. The 
scatter plots between cos#, cosines of the zenith angles 
of the charged leptons and fractional energies E^/E u of 
the charged leptons for diffrent directions of the incident 
neutrinos are given in Figs. 2 to 4. For a given Q 2 in 
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FIG. 2: The scatter plot between the fractional energies 
E^/Ev of the emitted muons and the cosine of zenith angles, 
cos#, for muon neutrinos with 1 GeV. The incident direction 
of the neutrinos is vertical. The sampling number is 1000. 
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FIG. 4: The scatter plot between E^/E^ and cos#. The 
incident direction of the neutrinos is diagonal. The other 
quantities are the same as in Fig. 2. The incident direction of 
the neutrinos is diagonal. The other quantities are the same 
as in Fig. 2. 
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FIG. 3: The scatter plot between E^/Ev and cos#. The 
incident direction of the neutrinos is horizontal. The other 
quantities are the same as in Fig. 2. 



Eq. (1), the energy Eg of the charged lepton and its 
scattering angle S is uniquely determined due to the two 
body kinematics. In Fig. 2, we give the case of vertically 
incident neutrinos {6 U = 0°). Here, as the zenith angles 
of the charged leptons are measured from the vertical 
direction in the SK case, the azimuthal angles of the 
charged leptons never influence over their zenith angle, 
and consequently the relation between their fractional 
energies and their zenith angles is uniquley determined 
as mentioned above. In Fig. 3, we give the case of 
horizontally incident neutrinos (6 U = 90°). Here, the 
azimuthal angle of the charged leptons has a potent 
influence on their zenith angle through the operation of 
Eq. (4). As is seen clearly from the figure, the cos# is 
widely distributed even to the backward for the same 
energy of the charged lepton. In Fig. 4, we give the 
intermediate case of the diagonal incidence (#„ = 43°). 

To connect our results with the analysis of the real 
experimental data , we finally need to take account of 
the energy spectrum of the incident neutrino in our 
calculation. For this purpose, we adopt the neutrino 
energy spectrum at Kamioka site obtained by Fiorentini 
et al.[6] and have carried out the following Monte Carlo 



procedure for a given cos Q v of the incident neutrino . 

Procedure A: we randomly sample the energy of the in- 
cident neutrino from the probability function which is 
composed of the combination of the neutrino energy spec- 
trum by Fiorentini it et al., which covers from 0.1 GeV to 
100 GeV at Kamioka site , with the corresponding total 
cross section for QEL. 

Procedure B: we decide whether the neutrino concerned 
should be attributed to the particle or the anti-particle by 
the random sampling from the corresponding probability 
functions given in Procedure A. 

Procedure C: we randomly sample Q 2 of the neutrino 
concerned for a given energy of the neutrino from Eq. 

(!)• 

Procedure D: we decide the energy of the charged lepton, 
Eg, from Eq. (3) and its scattering angle, 9 S , from Eq. 
(2). 

Procedure E: we randomly sample the azimuthal angle, 
(f>, by using a random number between (0,1). 
Procedure F: we decide the direction cosines of the 
charged leptons from Eq. (4), being accompanied by the 
scattering angle and the azimuthal angle obtained from 
Procedures D and E. 

We repeat a chain of the procedures A to F until 
we attain at the required trial number (ten thousands 
samplings per each case). 

In Figs. 5 to 7 for the three cases of 6 V (6 V = 0°,90° 
and 43°), we give the zenith angle distribution of the sum 
of /i+ and ii" for a given 9 U , obtained by Procedures A to 
F. If the SK assumption is valid, even if it may be of ap- 
proximation, the zenith angle distribution of the charged 
lepton should be of the delta function type with a peak 
around the direction of the incident neutrino. However, 
it is clear from the figures that the real zenith angle dis- 
tributions of the charged leptons deviate far from the 
delta function type distribution, and further the back- 
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FIG. 5: Zenith angle distribution of fi and fjT for and 
taking account of the overall neutrino spectrum at Kamioka 
site. The direction of the incident neutrino is vertical. 
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FIG. 6: Zenith angle distribution of /i + and pT for and v^, 
taking account of the overall neutrino spectrum at Kamioka 
site. The direction of the incident neutrino is horizontal. The 
other quantities are the same as in Fig. 5. 
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FIG. 7: Zenith angle distribution of /i + and fx~ for and , 
taking account of the overall neutrino spectrum at Kamioka 
site. The direction of the incident neutrinos is diagonal. The 
other quantities are the same as in Fig. 5. 



FIG. 8: Zenith angle distribution of e~ and e for u e and i7 e , 
taking account of the overall neutrino spectrum at Kamioka 
site. The direction of the incident neutrino is vertical. The 
sampling number is 10000. 
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FIG. 9: Zenith angle distribution of e~ and e + for u e and u e , 
taking account of the overall neutrino spectrum at Kamioka 
site. The direction of the incident neutrino is horizontal. The 
other quantities are thesame as in Fig. 8. 
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FIG. 10: Zenith angle distribution of e~ and e + for u e and z7 e , 
taking account of the overall neutrino spectrum at Kamioka 
site. The direction of the incident neutrino is diagonal. The 
other quantities are the same as in Fig. 8. 



ward scattering is too serious to neglect, which means 
rather large mixture of down going muon events into up- 
ward going events. It is also noticed that the azimuthal 
angles critically influence the translation from Fully Con- 
tained Events to Partially Contained Events (vicer versa) 
which is strongly dependent on their generation point in- 



side the detector. 

In Figs. 8 to 10, we give similar results for e~ and e + 
which correspond to Figs. 5 to 7. 

The charcteristics of the Figs. 8 to 10 for (anti- 
electron is essentially the same as in Figs. 5 to 7 for 
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(anti-)muon. 

V. CONCLUSION 

From these figures 5 to 10, it is surely concluded that 
the zenith angle distributions of the charged leptons un- 
der the SK assumption does not reflect the real zenith 
angle distribution of the incident neutrinos concerned. 
Of course, it is common sense that the direction of indi- 
vidual charged lepton cannot determine that of the neu- 
trino concerned uniquely, but only statistically. However, 
our present results clearly shows that the direction of the 
charged leptons under the SK assumption cannot dter- 



mine those of the incident neutrino, even if they accu- 
mulate experimental data statistically enough. Thus, it 
is surely concluded that we could not extract any defi- 
nite conclusion around the neutrino oscillation from the 
measurement on the charged leptons under SK assmp- 
tion which occur inside the detector and is regarded as 
the most reliable one from the point of the experimental 
accuracy. 

The relation between the zenith angle distributions of 
the charged leptons which are really measured experi- 
mentally and those of the parent neutrinos, taking ac- 
count of the scattering angle and the azimuthal angle of 
the charged leptons correctly, will be shown in a subse- 
quent paper. 
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